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We develop a general approach to calculating "nonuniversal" prefactors in static and dynamic 
correlation functions of ID quantum liquids at zero temperature, by relating them to the finite 
size scaling of certain matrix elements (form factors). This represents a new, powerful tool for 
extracting data valid in the thermodynamic limit from finite-size effects. As the main application, 
we consider weakly interacting spinless fermions with an arbitrary pair interaction potential, for 
which we perturbatively calculate certain prefactors in static and dynamic correlation functions. We 
also non-perturbatively evaluate prefactors of the long-distance behavior of correlation functions for 
the exactly solvable Lieb-Liniger model of ID bosons. 
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I. INTRODUCTION 

One-dimensional (ID) quantum liquids of bosons, 
fermions and spins are conventionally described using an 
effective hydrodynamic approach known as the Luttinger 
liquid theory ^[^. This theory predicts the long-range 
behavior of equal-time correlation functions at zero tem- 
perature, which one obtains as a series expansion with 
power laws controlled by a dimensionless Luttinger liquid 
parameter K > 0, see Eqs. ([I])- ([3]). While the "univer- 
sal" parameter K is related to thermodynamic properties 
and can be easily extracted from numerical or exact solu- 
tions, the "nonuniversal" prefactors in the series expan- 
sion, e.g. Am, Bm,Cm, sec Eqs. dB-dSt are usually not 
known except for a few cases [a-flOj. At the same time, 
these prefactors set the actual scale of observable corre- 
lations, consequently determining them is an important 
theoretical challenge. 

In this article, we first develop a general technique 
for calculating these nonuniversal prefactors by combin- 
ing the Luttinger liquid Hamiltonian with the analysis 
of the finite-size properties of certain matrix elements 
(form factors). We then consider dynamic response func- 
tions such as the density structure factor and the spectral 
function, see Fig. [TJ It has been shown recently |ll| - [l9| 
that dynamic response functions generically have singu- 
larities which can be described by effective Hamiltoni- 
ans of impurities moving in Luttinger liquids. Analysis 
of the finite-size properties of these effective Hamilto- 
nians allows us to extend the approach to various dy- 
namic response functions. To demonstrate it, we per- 
turbatively evaluate several prefactors of equal-time cor- 
relation functions and dynamic response functions for 
weakly-interacting fermions. We also calculate, non- 
perturbatively, various prefactors for the exactly solvable 
Lieb-Liniger model [lo, [HI of ID bosons with contact 
interactions. The latter model has been realized with ul- 
tracold atomic gases [l^, and its correlation functions 
can be measured using interference f23l |2^ . analysis of 
particle losses [1^, photoassociation [25|, or Bragg and 



photoemission spectroscopy [26|. 

This article is organized as follows. In Sec. II A, we use 
linear Luttinger liquid theory to work out the connec- 
tion between prefactors of equal-time correlation func- 
tions and lowest energy form factors. In Sec. IIB, we 
show that the relative spectral weights of all low energy 
form factors can be fixed based on universal nonlinear 
Luttinger liquid theory [16]. In Sec IIC, we show that 
the effective field theory of impurities moving in Lut- 
tinger liquids allows to extend the relations between form 
factors and prefactors to dynamic response functions. In 
Sec. Ill we present the results of the perturbative calcula- 
tions of various prefactors for weakly interacting spinless 
fermions. We summarize our results in Sec. IV. Some of 
the technical details are contained in the Appendices. 



II. RESULTS FROM EFFECTIVE FIELD 
THEORY 

A. Prefactors of equal-time correlators from the 
Luttinger liquid theory 



The Luttinger liquid theory [l|-l5| predicts the behav- 
ior of the correlation functions for spinless bosons and 
fermions of density po at pqx ^ 1 as (here kp = npo) 
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Here p is the density operator, and ipp{ipB) is the 
fermionic (bosonic) annihilation operator. The Hamilto- 
nian describing these correlations is written as (we follow 
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notations of Ref. [1]) 



(4) 



where v is the sound velocity, the canonically conju- 
gate fields (/){x),9{x) have the commutation relation 
[(j){x),V9{x')] = iTTS{x — x'), and the components of the 
fermionic (bosonic) fields with momenta (2m + 1/2 ± 
\/2)kF are written as 



1pF(B){x,t) 



J{2m+l/2±l/2)[kFX-rl>{x,t)]+ie(x,t) 



(5) 



For repulsive bosons, one has K > 1, while for repulsive 
(attractive) fermions K < 1(> 1). 

One of the reasons for the success of the Luttinger 
liquid theory is its ability to predict certain finite-size ef- 
fects [3, [131 due to the conformal invariance of the Hamil- 
tonian Below we will show that conformal invariance 
can also be used to relate nonuniversal prefactors to the 
scaling of certain form factors [see Eqs. (ITT]) - p^ . (^5]) ]. 
which constitutes the main result of this article. Form 
factors can be evaluated perturbatively or numerically 
for finite-size systems, and are known for certain inte- 
grable models such as the XXZ [1^, [29Ljthe Calogero- 
Sutherland [sO], and the Lieb-Liniger 2ll, [sij models. 
Thus relations PT|) - ([T^ . (P5|) provide a powerful tool with 
which one may interpret finite-size effects and make pre- 
dictions which are valid in the thermodynamic limit. 

Let us start by considering interacting spinless 
fermions. Using the resolution of the identity in the ex- 
pectation value {ipp{x,t)ipp{0)), we get 



{^pUx,t)4'Fm = Y.<'^'^'''~ 



{k,u\i;F\N) 



(6) 



where {k,uj\ipp\N) is a form factor of the annihilation 
operator, \k,u!) denotes an eigenstate of — 1 particles 
with momentum k and energy cj, and |A^) is the ground 
state of N particles. For simplicity we assume N is odd 
so the ground state is non-degenerate. For a finite sys- 
tem, k and uj are not continuous, but will be quantized 
and consequently the spectral function is a collection of 
delta functions in (fc,cj). We will now obtain a similar 
representation from the Luttinger liquid theory and com- 
pare it with Eq. ([6]) to obtain the nonuniversal prefactors 
Cm- Hamiltonian (j4]) can be written using left- and right- 
moving components ^Pl{r) = ± (p/\/K, Q which 
dictates the time dependence of the 
nent of (V3|.(a;, t)V'F(0, 0))//9o at pQ\x±vt\ ^ 1 as 



compo- 



pZ(2m+l)fcFa; 
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2z(-l)" {i{vt + x) + Of^'^ {i{vt -x)+ Qf^' 
where 

PfmR) = (2"i + l)^^s:/4 ± (2m -f l)/2 + \/AK > 0. (8) 



The coefficients Cm appeared in Eq. ([7]) due to the com- 
parison of the t = limit of {'ipp{x,t)ipF{0,0)) / Po with 
the right hand side of Eq. ([3]). The two factors in the 
denominator of Eq. Q describe contributions from left 
(right)-going excitations near the Fermi points which 
propagate with velocities ^v. The signs of the infinites- 
imal shifts in the denominators ensure that only excita- 
tions with negative (positive) momenta can be created at 
the respective branches. For a finite system with periodic 
boundary conditions on a circle of length L, conformal 
invariance dictates (see e.g. Ref. [3|) that Eq. (O gets 
modified as 
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We can now expand the terms in the parentheses in a 
Fourier series 
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Comparing Eqs. (Pjl- dTUl) to the right hand side of Eq. (0), 
one can clearly identify contributions from low-energy 
and momenta particle-hole excitations at the right (left) 
Fermi branches with energies 2nvn±/L > and mo- 
menta ±2nn±/ L. Additionally, m inter-branch pairs of 
momentum 2kp each are formed by successively deplet- 
ing discrete states below the left Fermi point (m > 0) 
and occupying the lowest-available states above the right 
one. On top of that, an additional hole is formed on the 
left branch, giving a total contribution of {2m + l)kp > 
to the momentum. The contribution from n+ = n_ — 
gives the scaling of the "parent" form factor 
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where |m, — 1) is the lowest energy state of — 1 
fermions with momentum (2m -|- 1)^^;-. For Galilean in- 
variant systems, states with different m can be obtained 
from the ground state by a center of mass Galilean boost. 
We see that as a consequence of the criticality of the Lut- 
tinger liquid, form factors of the annihilation operator 
have nontrivial scaling with the system size, and the pref- 
actors in front of these nontrivial powers of L are directly 
related to the prefactors in the correlation functions. For 
density correlations and bosons, similar relations can be 
worked out, and are given by 
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FIG. 1: (Color online) (a) Spectral function Aik^iS), with 
shaded areas indicating the regions where it is nonzero, and 
notations for prefactors. (b) "Parent" state responsible for the 
singularity at « — £(fc) > in spectral function: it contains 
a hole corresponding to a mobile "impurity" at fcj, « — fc, 
and one particle at each Fermi point. We also illustrate the 
finite size quantization of the momenta of the impurity and 
of excitations at the right Fermi point. 



Eqs. ([TT|) - ([T3l) allow one to evaluate the prefactors in 
Eqs. (HJ-® by identifying a single, simplest "parent" 
form- factor for each of the operators p,i)B and iI)f, re- 
spectively. Results for bosons are simply generalized to 
describe Luttinger liquids of spins on a lattice with stan- 
dard substitutions [4] , and in particular Eq. ()13|) explains 
the coincidence noticed in Refs. for the spin- 1/2 

XXZ model. 



B. Distribution of spectral weights among low 
energy form factors and the multiplet summation 
rule 

Let us now consider the spectral weights at finite en- 
ergies and momenta in the vicinity of the Fermi points 
(for concreteness we will discuss fermions). In a finite- 
size system, the spectral function is a collection of delta 
functions in (fc, w), weighted by form factors, see Eq. ([6]). 
From field theory, we not only determine the parent 
form factor, but also form factors associated with states 
containing low-energy excitations on top of the parent 
state, described by nonzero n+ and n_ .Within the lin- 
ear spectrum approximation near the Fermi points, for 
n± > 1 certain particle-hole states are degenerate, thus 
within oc 1/i accuracy, from Eq. ^ and Eqs. (P)- p^ 
the total spectral weight which falls into the degener- 
ate subspace ( "multiplet" ) with given , n_ can be ob- 
tained from the parent form factor by multiplying it by 
C{n^, ^F^ii)C{n-, ^p^l). However, within cx accu- 



racy, states with sufficiently large n+ and n_ are not de- 
generate due to the nonlinearity of the generic fermionic 
spectrum, as can be illustrated by the case of weakly in- 
teracting fermions. Then C{n^, fipji)C{n-, fip^^) gives 
only the total spectral weight within each multiplet, and 
the conventional linear Luttinger liquid theory doesn't 
distinguish the oc 1/L^ splitting of the contributions 
within each multiplet. However, the splitting of the spec- 
tral weights within each multiplet is also universal, and 
can be understood based on the universal theory of non- 
linear Luttinger liquids developed recently [la]. Below 
we will illustrate such splitting for the case of n_ = 0, 
i.e. when only excitations at the right branch are created. 

Within the nonlinear Luttinger liquid theory, we evalu- 
ate various correlators (e.g. Eq. (IT]),([3])) by first express- 
ing the fermionic creation and annihilation operators in 
terms of fermionic quasiparticle operators V'_r(l)- The 
fermionic operators are related to the quasi particle op- 
erators using 

Taking the nonlinearity in the spectrum of the original 
fermions into account [l6j . we obtain the Hamiltonian 
Hii + Hl in terms of the quasiparticle operators with 
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where m* is the effective mass, characterizing the nonlin- 
earity of e(fc). Its inverse can be expressed through low 
energy parameters as [32j 
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with h denoting chemical potential and v the sound veloc- 
ity. Expressions for correlators will include contributions 
from the left and right Fermi point with given /iF,fl(L)i 
i.e. terms 



oc \^exp 
exp 



2m^^p,ji I dy : ^\{y,t)ipR{y,t) : 
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where the average is taken over the infinite chiral Fermi 
sea |FS) with all negative momenta occupied, and nor- 
mal ordering is with respect to this vacuum. A similar 
contribution from the left-movers will also appear. 

In a finite system, momenta of the fermionic quasi- 
particles are quantized near the Fermi point in incre- 
ments of 2n/L. We denote by pi > P2 > ■■■■ > Pn > 
particle excitations carrying momenta 2Trpi/L, and by 
qi < (72 < ••■ < Qn < hole excitations carrying mo- 
menta —2TTqi/L. Then 2Tm+/L ~ 27: / L'Y^^pi — qi is 
the total momentum of particle-hole excitations near the 
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right Fermi point. Within the linear spectrum approxi- 
mation, all states with the same n_|_ are degenerate, while 
due to nonlinearity each of them acquires an energy shift 

= J^E(p'+9')- (18) 

i 

We can introduce a complete set of intermediate states 
in Eq. (|17p between the two exponents and organize them 
by the momenta of particle-hole excitations as follows 
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where in the first line we have moved the posi- 
tion and time dependence of the operators over to 
the states, and in the second line we have Fourier 
transformed the creation and annihilation operators 
and performed the integral over y. The correlator 
above can be exactly evaluated using the methods of 
Ref. [33]- There they consider a "boundary state" 

exp -(a + TO)X;fc^o (^EpV'jiV'p+fc) |0), which we 
can identify as the one obtained by action of the string 
operator in Eq. (|19p on the infinite chiral vacuum if we 
map m + a ^ y/T^FM- Consequently we obtain the fol- 
lowing rcsuh from Eqs.(59)-(61) of Ref. [H]: 

f{{Pt,Qi}) 
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r(i + V7il^)r(-g) 

Normalization of the spectral weight leads to the fol- 
lowing "multiplet summation rule" (see Appendix A for 
an explicit demonstration) 

E \f{{p^.q^})? ^C{n+,^lp^n)■ (22) 

Contributions from the left Fermi point are accounted 
for similarly, and the total form factor is a product of 
these two terms. 



C. Prefactors of dynamic response functions from 
effective theory of impurities moving in Luttinger 
liquids 

We now consider the dynamic response functions: the 
density structure factor 

S[k,oj)^ [ dxdf^'^'-''^'>{p{x,t)p{0,0)), (23) 



and spectral function A{k,uj) — —^lmG{k,uj)signuj 
where the Green's function G{k, ui) is defined as [3J| 

Gik,uj) = -i I da;dte*('^*-'=^)(r[v3(x,i)Vi^(0,0)]). (24) 



To be specific let us consider the spectral function for 
w > and —kp < k < kp- In addition to the Fermi 
points, the field theoretical description of the singular- 
ity at —£{k) > (see Ref. 17[; we follow the notations 
contained therein) involves a mobile "impurity" with mo- 
mentum kh w — fc moving with velocity Vd — de{kh) / dkh, 
see Fig. [T] For non-interacting fermions, any spectral 
weight is absent at a; > and —kp < k < kp, and for 
weakly interacting fermions the configuration responsi- 
ble for a feature at lu « —£{k) is illustrated in Fig. [1] it 
has one particle at each Fermi points, and a hole corre- 
sponding to the impurity. While for stronger interactions 
such a simple interpretation of the impurity is absent, the 
field theoretical description still remains valid The 
Hamiltonian used in this approach takes the form 



Hd = J dxd^a 



e[k) - ivd^ 
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Here the operator d{x) creates a mobile hole with mo- 
mentum k and velocity Vd = de{k)ldk. The interaction 
Hamiltonian describes the impurity interacting with the 
left and right movers of the Luttinger liquid. 

The spectral function A{k,uj) in the vicinity of — e(fc) 
is written as 



A{k,u:) cx J dxdte'^'{d^x,t)d{0,0))H^^+H,+H,„, 
= Ao^+{k) I dxdte'^"*D{x,t)L{x,t)R{x,t), 



(26) 

where Suj ~ u + s{k), D{x, t) — 5{x — Vdt) is the impurity 
correlator, L{R){x,t) = {i{vt±x)+Oy°-+'^'^^ [SSj, and 
we introduced a prefactor Ao^^{k) which will be deter- 
mined by a comparison to the form factors. After the x, t 
integration, Eq. results in 



A{k,uj) = e{5uj] 



2TTAo.+{k)5uj-f"'- 



r(l - Pq,+){v -t- Vd)^'°^+•^ {v - VdY»-+'^ ' 
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In finite-size systems, L{x, t) and R{x, t) get modified, see 
Eq. (dni). The change of D{x,t) to J2no e^"""^^-'''^*)/^ 
corresponds to the quantization of the impurity momen- 
tum. At fixed k, the shift of the momentum of the impu- 
rity can be expressed a.s no = n_ — n+. Combining these 
terms, we get 



A{k,uj) 



^ siSLj-AE-?^{v-Vd) 
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L 



iv + Vd) X 



(2^) 



2 — ^0,4 



-C(n+,^o,+,i?jC'(n_,^o,+,L), (27) 



where AE is a universal oc 1/i shift of the edge posi- 
tion [l^. Thus the finite size structure of the response 
function is given by the sum of two generically incommen- 
surate frequency "ladders" at arbitrary k, in contrast to 
the vicinities of Fermi points. Analysis of the scaling of 
the parent form factor with = n_ = then leads to 



{k;N + l\i^UN) 



2— M0.+ 



(28) 



where |fc; + 1) denotes the lowest energy state of A^ -|- 
1 fermions with total (discrete) momentum k. Similar 
relations can be derived for the density structure factor 
and the boson spectral function, and for each fc, the left 
hand side consists of a single form factor which connects 
the ground state with the lowest energy state at total 
momentum k, while the right hand side shows scaling 
with the exponents of Ref. [17]: 



5*0 (fc) /27r 
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We note that in Eqs. dSH!)-® k has to be fixed before 
taking the limit L — >■ cxi, since e.g. the k ^ kp and 
L — ^ cx) limits do not commute as has been shown in 
nonlinear Luttinger liquid theory [l3, EB, [l3l ■ 

Eqs. (ITT])-([ni),(I2S])-(I2ni) rely on the structure of the 
low-energy excitations at a given momentum, prescribed 
by the field theory; they are valid for all Luttinger liquids 
in ID irrespective of microscopic interactions and can 
be used as a convenient tool to interpret the results of 
numerical studies. Below we illustrate their power by 
obtaining new nontrivial results for weakly interacting 
fermions and also present numerical data on some exact 
results |40| we obtained from the analysis of the finite 
size form factors of the Lieb-Liniger model. 



III. PERTURBATIVE CALCULATION OF 
FORM FACTORS FOR WEAKLY INTERACTING 
SPINLESS FERMIONS 

We use the non interacting Fermi gas in ID as our 
unperturbed state. The ground state for a system of A^ 
non interacting fermions of mass M occupying a length 
L with density pQ — N/L is characterized in momentum 
space by N real momenta {—kp, ...,kp}, increasing by 
increments of 2tt/L, with kp = tt{N — 1)/L. Here we 
have assumed that A^ is odd so that the ground state is 
non-degenerate. 

Defining the ground state |FS) to have resulted from 
the action of A'^ creation operators on an empty vacuum 
starting from the left-most momentum gives us a conven- 
tion to specify the relative phases of various states that 
will be used in the calculations to follow. Moreover in 
the subsequent calculations it is only the relative phases 
between states that is pertinent since we are interested 
in absolute squared values of the form factors. 

To this system we add a weak four fermion interaction 



where V{q) is the Fourier transform of the pair interac- 
tion potential V{r). 

For weakly interacting fermions we may directly eval- 
uate form factors in the left hand sides of Eqs. (fTT|) . (|13p . 
and ([28|) using conventional perturbation theory and 
extract prefactors. For example, since iio,+ = — 1 + 
O(y^) [m, one can expand the right hand side of 
Eq. (pS)) in powers of V as 



(27r)3Ao,+ (fc) 
L3 



X [1 + (^0,+ + 1) log(L/27r) + 



(32) 



While in an infinite size system this expansion is not 
convergent, for finite L it is well defined if one keeps L 
finite but large, and then takes the limit V{r) — >• 0. We 
treat the interaction term perturbatively, and can e.g. 
write the expansion of the ground state as 



|Ar) = |FS)+^ 



|a) 



(«I^|FS) ^ 

Eps — Ea 



(33) 



where |FS) denotes a filled Fermi sea. Similar pertur- 
bative expressions can be written for the states on the 
left hand sides of Eqs. (im) . ([T5)) . and (1^51) . and one can 
then straightforwardly evaluate the scaling of form fac- 
tors with integer powers of L. Due to momentum con- 
straints, only few intermediate states contribute within 
lowest order perturbation theory. Eg. for AQ^+{k) (see 
Fig. [T]a), the only sequence of states which contributes 
is the following: first, V creates two particles at the right 
and left Fermi points, and two holes at -t-fc and — fc; (see 
Fig. H^), second the operator V'^(O) — -^'^pi'Hp) 
in a hole at +fc, and we end up with the final state 
\k.N + l). 



6 






2fcf + k 



A. Calculation of Ao,+ {k) 



We expect the spectral function in the region co > 
0,k £ (— fc_F)(see Fig. [1^) to behave as 
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FIG. 2: The states (a) - (c) shown above are the only inter- 
mediate states that give non-zero contributions in the pertur- 
bative calculation of the form factors used to determine the 
prefactors Ao.+ {k), Ai^-{k) and Ai^+{k), respectively. 



where As captures the shift in the dispersion curve of 
a non interacting fermionic system when interactions are 
introduced. To first order the shift is obtained by consid- 
ering the Hartree-Fock corrections to the energy of eigen- 
states. The exponents fio^+, fJ.o,+,L, fJ-o,+,R (see Eq. (PT|) 
and Ref. (sH]) to leading order can be obtained using: 
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\ nkp 7r(fc — kp) Tr{k + kp) J 

l + OiV). 
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(35) 
(36) 



Using Eq. (|28|) and the value of /io,+ in the limit of 
vanishing interactions from Eq. p6p . we may obtain the 
prefactor ^o,+ (^) from: 



acting system perturbatively in V in the following way: 



(87r3) 



|(fc,iV-M|V^(0)|iV)|2. (37) 



|iV) 

|fc,A^ + l) 



|FS)+El^l")+0(n 



Eps — Ea 



Thus we need to evaluate the form factor 

(fc,7V + l|^t(o)|7V) = 1 (fc,7V + i|^^t^,|iV), (38) 
V p„ 

where |fc, + 1) is the lowest energy eigenstate of the 
interacting system with N + 1 particles and total mo- 
mentum k and \N) is the N particle ground state. Since 
there is no exact answer for such a form factor for a 
generic interacting Fermi gas, we expand the ket vectors 
in Eq. ([55)) in terms of the eigenstates of the non inter- 



la) 



El., 



k,N+l 



En 



0{V^), 



(39) 



where, the state |FS) is a filled Fermi sea of N particles, 
and the state \k,N + 1)'°) is an eigenstate of the N + 
1 non interacting fermions, that has two particles with 
momenta —kp — 2tt / L and kp + 2n/L, and a hole of 
momentum ~k, on top of |FS). From now on we will 
specify various states of the non-interacting system by 
describing the configuration of the state with respect to 
|FS). The sum over \a) runs over the entire Hilbert space 
of the non interacting system except the zeroth order 



7 



state on the right, and Ea is the energy of the eigenstate 
I a) of the non interacting system. 

When we use the expansions in Eq. (j39]) to hnear order 
in expression psp for the form factor, we obtain no con- 
tribution from terms that are zeroth order in V since it is 
not possible to create the state \k,N + 1) by acting on 
|FS) with just a single creation operator. When we con- 
sider the two terms generated by taking the zeroth order 
term from one of the ket vector expansions and the first 
order term from the other it becomes possible to obtain 
non zero contributions. Thus we are left considering two 
terms: 

x{a\ ^ V{q)^pl^^4>l, _^Mp\FS} (40) 

Q-P-.P' 

and 

x(fc,iV + l|W J2 V{q)i>l+^iJ:l,_^i>p4p\a)m 
q,p,p' 

I 



We see that only the first term leaves a nonzero con- 
tribution for the following reasons. In term (j40|) . the 
matrix element of the momentum conserving interaction 
V limits the possibilities for la) to be only states with 
zero momentum. On the other hand, for the matrix ele- 
ment of the creation operator to be non zero we require 
that the state \a) must have the same configuration as 
|fc,iV -I- 1)(°) but with one more hole. Since the state 
|fc, A^-|- 1)(°' already has a hole, the only admissible state 
is one with two particles of momenta kp + 27r/L and 
—kp + 21: / L and two holes of momenta k and —k and is 
shown in Fig. [Ifa). When we consider term (HIT) we see 
that the states \a) which will give non zero contributions 
must have total momenta k and must contain one ad- 
ditional particle over the ground state. There can exist 
no such state since |fc| < kp, and consequently we may 
disregard term ((4T|) . 



Thus the only contributions we are left with are due 
to the term in expression for the state in Fig. [UJa): 



{k,N + l\ij\Q)\N) 



[kp + 2n/Lf e 

M M ' 

M{V{kp -k + 2tt/L) - V{kp + k + 2n/L) - V{-kp -k- 2tt / L) + V{-kp + k- 2tt / L)) 

2L^I^[{kp + 2TT/LY -k^] 
+0(y2). (42) 



Using the correspondence between the prefactor 
Ao,+ (fc) and the form factor given in Eq. p7p we obtain 



^o,+ (fc) 



L3 



B^3) 



\{k,N + l\i,\Q)\N)\'' 



M'^{V{kp + k) ~ V{kp - k) f 



/c2)2 



(43) 



to leading order in V and where we have used the fact 
that V{r) is real and symmetric, thus V{—q) — V{q). 
The divergence as k — >■ ±kp from the denominator of 
Eq. still leads to a finite integral over ui > when 
we substitute the expression for Ao,+ in Eq. (j34p . 



B. Calculation of Ai_(fc) 

In the region w < 0,fc S {kp,3kp) (see Fig[T^) the 
spectral function behaves as 



A{k,Lo) w 27rAi,_(fc) 



kip 



{k~2kpf 



2m 

2m 



Ae 
K Ae 



r(i - iii,-){v + vdT'--^ {v - vdT'---'' ' 

{V{k' -k + 2kp) -V{k' -kp)) 



27r 



+0{V^). 



(44) 



The exponents /ii___L, Mi.-.j?, are given by Eq. ((3T|) 
and Ref. (ssj . and can be written to leading order as 
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= -1 



K- 



2tt 



m{V{0) -V{2kF)) m{V{kF + k) -V{0)) 
inkp 2'K{k + kp) 

m{V{0) -V{2kF)) m{V{kF -k) -V{0)) 



47rfcF 2-K{k-kp) 
m(y(0) - F(2A:f)) m(F(fc - fej.) - F(0)) m(F(fc + A:f) - V(0)) 



Trfci 



7r(fc — kp) 



■K(k + fcj 



P-i,- = -i + oiv). 



(45) 
(46) 



r 



Using the value of the exponent in the absence of in- 
teractions and the finite size version of Eq. (j44p we may 
write 



-\{k,N-l\m\N)? 



Consequently we need the matrix element 



(47) 



{k,N-imO)\N) = 



{2kF 

1 

7l 



kh,N~l\ij{0)\N) 



{2kF + kh\J2^Jp"\N), 



(48) 



where kh ^ k — 2kF G {—kF,kF), and the state \2kF + 
kh, N — 1) is the lowest energy eigenstate of momentum 
2kF + kh oi N — 1 interacting fermions and \N) is the N 
particle ground state of the system. Again we perform 
a perturbative expansion of the ket vectors in orders of 
V{q) in eigenstates of the free fermions as in Eq. (p9| . We 
note that the unperturbed version of the excited state, 
\2kF + kh, N — will contain two holes with momenta 
kF,kh and a particle of momentum —kp. 

We find that there is no zeroth order contribution since 
the state \2kF + kh, N — l)^*^^ cannot be obtained by the 
action of a single annihilation operator on the ground 
state, since it contains two holes and a particle over the 
ground state. Thus we look to the first order terms of the 



form given in expressions PO)) and (|^T|) . In order for the 
linear order terms to give a non-zero contribution we need 
the matrix elements of V , which can only connect states 
of the same momentum, and that of the annihilation op- 
erator between (a| and |FS), to simultaneously be non- 
zero. The first requirement automatically narrows down 
our choices to two subsets of intermediate states {|ck)}. 
One subset contains states of zero momentum while the 
other contains states of momentum 2kF + kh- For the ma- 
trix element of the annihilation operator to be non-zero, 
the states \a) of momentum 2kF + kh must be identi- 
cal to the filled Fermi sea but with a single additional 
hole - however this is impossible since the momentum of 
such a state can at most be fc^?- Thus we must focus our 
attention on the term 



E 

|a) 



2L3/2(i5ps-i?„) 



{2kF + kh,N ~l\Y,i'p"\a) 



p" 



{a\ ^ Viq)^jl+Jl,_gi^p4p\FS). 



The only state |a) that gives a contribution to this term 
contains two holes of momenta kp and kh and two parti- 
cles of momenta —kp — 2Tr/ L and 2kp + kh + 2Tr/ L, and 
is shown in Fig. [l{b). Consequently the matrix element 
to first order in V can be calculated as follows: 



We may then express the prefactor v4i^_(fc) in terms Eq. P5|) : 
of the matrix element using a similar correspondence to 

AiMk) = j^\{k,N-l\m\N)\' 

^ M^(V{2kp)-V{kp-k)f 

2,2TT^kl{kp - kf ' ^ ' 
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where we have substituted kh = k — 2kp to obtain the 
final answer in terms oi kp < k < 3kp. 

As in the case of ylo.+ (fc), we note again that the diver- 
gence as fc — >■ kp from the denominator of Eq. (15(11) still 
leads to a finite value when we substitute the expression 
for Ai^-{k) in Eq. (|44l) and integrate over u < 0. 



Moreover we obtain perturbative expressions for the 
exponents fii^^^n, ^i^+^l and from Eq. (|3T|) and 

Ref. [H: ' ' 



C. Calculation of Ai^+{k) 

When oj > 0, fc e {kp,3kp) (see Fig. [1^) the spectral 
function behaves as 



A{k,uj) « 27rAi,+ (fc) 



kj,-{k- 2kpf 



2m 



Ae 



Ae 



Ae = 



■ r(l - ni.+){v + VdY^-+-^ {v ~ w<j)^i.+.« ' 
• Ay 

— {y{k' -k + 2kp) - v{k' - kp)) 

27r 



(51) 



Ml, + ,L 



"I MkTkp) ^ ^ 



2m(y(fc + fcj' 



m(T/(0) - F(2fcF)) , m(y(fcF - fc) - V(0)) , ^^^y^. 



n{k + kp) 



Auk I 



0{V^). 



27r(fc - kp) 



= -3 + 0(T>). 



(52) 
(53) 



Using the value of the exponent in the absence of in- 
teractions and the finite size version of Eq. (j5ip we may 
write 



L5 



-\{k,N -l\i^\Q)\N)\ 



(327r5)' 

Thus we first focus on the matrix element 



(54) 



(fc,A^ + l|^^(0)|A^) = {2kp~kh,N + 1\'^\Q)\N) 
= l={2kp~kn\Y,i^l,.\N), 



(55) 

where kh = 2kp — k € (—kp,kp), and the state \2kp — 
kfi,N + 1) is the lowest energy eigenstate of the N + 1 
particle interacting system with momentum 2kp — k}i and 
\N) is the N particle ground state. We again expand the 



ket vectors in Eq. ([55]) as in Eq. ([M]) . The corresponding 
eigenstate of the free fermions is \2kp — kh,N + 
an TV -|- 1 particle state with two particles with momenta 
kp + 2tt / L and kp + An/ L, and one hole with momentum 
kh ■ Here too we must perturbatively expand in the states 
for which we calculate the above matrix element. There 
is no zeroth order contribution since the state \2kp — 
kh,N + cannot be created by the action of a single 
creation operator on the ground state. 

Upon considering the first order terms of the form 
given in expressions (|40p . (|4ip we find the only non-zero 
contribution is to the term 



E 



1 



[a 



2Ly^E2k,-k,,N+i - E^) ^„ 

X {2kp-kh,N + l\ J2 ViQ)^^l+g^^l'-Jp4p\a) 
q,p,p' 
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from the state depicted in Fig.[2l^c). Due to the kinematic 
constraint imposed by V, the states that may contribute 
to the matrix element must either have zero momentum 
and must be connected to the state \2kp — kh, by 
the action of the creation operator, or have momentum 
2k p — kfi and must be connected to the ground state 



by the action of the annihilation operator. The latter 
condition is the only kinematically feasible one, and the 
only allowed state is the one in Fig. [SJc). 

Thus we may collect all the terms that are linear order 
in V as follows: 



{kp+4n/Lf , {kp + 2n/L)^ kl i2kp - k^ + 6n/L)^ _ {kp - k^f , r.ni r , 
E2k.-,,.N+i-E,, - — + — — — +0(1/L), 

/9i. y Ma^MlUnMM\ M{V{kp - kn + - V{kp - kn + 27r/£)) 2^MV'{kp - kn) 
{2kp-k,.N^l\^HmN) - L3/2(fc^ _ fc,)2 L^I\kp-k,Y ■ 



Note that in the above expression, the appearance of the 
derivative is essential to give the correct power of 1/L, 
which is evident when one considers the correspondence 
between the prefactor Ai,+ (fc) and the form factor ([55)1 
given below. 

From finite size scaling the prefactor +(fc) can be 
obtained using the matrix element above: 



327r' 



■|(fc,iV + l|V-t(0)|7V)| 



2 _ APVik-kpY 



%TT^{k-kpY ' 

where we have substituted k^ = 2kp — k to obtain the 
final answer in terms of kp < k < 3k p. We note that 
the divergence of Eq. ([571) as k ^ kp is generically the 
same oc (fc — kp)~'^ divergence seen in Ai^_(fc — >■ kp). 
Yet again the integral of the spectral weight over uj > 
remains finite as can be seen from Eq. (j5ip . 

The procedure to obtain the prefactors An>i,±{k) is 
similar to the one used so far. In general, to lowest non- 
vanishing order A„>i^±(fc) will be cx V^^", but one needs 
to sum contributions to the form factor expansion from 
a rapidly growing number of intermediate states. More- 
over terms appearing with higher orders of V need to 
be carefully separated from terms cx (log(L))" that gen- 
erate corrections to the exponent (see Eq. (15^ ). The 
latter procedure is demonstrated in Sec. IIIF where we 
calculated Ao._ (fc) to the lowest non- vanishing order in 
V. 



D. Calculation of prefactor Ci 

We now consider the prefactor Ci defined in Eq. ([3|). 
We denote by |m = 1,A^— 1) the lowest energy A^ — 1 par- 
ticle eigenstate of the interacting fermion system of mo- 
mentum 3kp. The corresponding state of the free fermion 
system (zeroth order term of the perturbative expansion) 
has two holes at the left Fermi point and one particle at 
the right Fermi point and is depicted in Fig. ^a). We 
wish to calculate the annihilation operator form factor 





FIG. 3: The state in (a) depicts the parent state with respect 
to which we calculate the annihilation operator form factor 
to obtain prefactor Ci. That state in (b) depicts the only 
intermediate state which gives non-zero contribution to the 
form factor in the perturbative expansion to first order in V . 



for this state and the N particle ground state \N): 
{m^l,N~l\i;{0)\N) = ^(m = l,iV- 1| ^ V'p"!^)- 

(58) 

We may again expand the states above like Eq. p9l) . 
We find that the first non-zero contributions in the per- 
turbative expansion appear at the first order in V. To 
this order we obtain the following two terms by expand- 
ing the states in the bra- and the ket vector perturba- 
tively in V: 



-j,,5iji--^(„..i,«-iii:.v 

\a) p" 

and 

1 



(59) 



E 



\a) ^ ■ p" 



al^VlFS) 



= \,N-\\Y, V{-q)^l+^^l,_^Mp\<^)- 



(60) 
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The states which give non-zero contributions to the 
term (|59p must have momentum. Furthermore we must 
be able to transform them into the state \m — 1,N—1) by 
the action of a single annihilation operator. This leaves 
only one possibility shown in Fig. [3jb). While for the 
term (|60p to be non-zero we require states which have 



momentum 3kp which can be obtained from the ground 
state by the action of a single annihilation operator, i.e. 
a state with one hole of momentum —Skp, which is im- 
possible to achieve. 

Thus we find the form factor as follows: 



4 , (fcF - 2^/i)2 (kp + 27r/L)2 (3fcf )2 4kF (, , tt 

i^FS — -C'^h = — n^P H 

2M 2M 2M 2M M \ L 

I 1 AT il/rmipq\ M{V{2kp + 21, /L) - V{2kF)) M7rV'{2kr) 

{n,^l,N-l\m\FS) = = ___. (61) 



We may express Ci in terms of this form factor using 
the correspondence we have derived between the form 
factor calculated above and our desired prefactor: 

(62) 

E. Calculation of prefactor A2 

To obtain the prefactor A2 in Eq. ([T]) we must calculate 
the following matrix element: 

(m = 2, N\p{0)\N) ^{m^2,N\^Yl ^l'+q'^k'\N), 

q',k' 

(63) 

where jm = 2, N) is the lowest energy eigenstate of the in- 
teracting system with momentum Akp. The unperturbed 
version of this state, for TV free fermions, contains two ad- 
jacent holes in the vicinity of the left Fermi point and two 
adjacent particles in the vicinity of the right Fermi point. 

We evaluate the form factor in Eq. (|63|) by again per- 
turbatively expanding the ket vectors in the free fermion 
basis, in powers of V as in Eq. (l39l) . We again disregard 
the terms that are zeroth order in V since the density 
operator cannot create the m = 2 Umklapp state from 
the ground state, and consider the first order terms for 
the first non vanishing contributions. 

The first term of 0{V) is 

|a> q',k' 

{a\ J2 V{q)^l+,^l,_^i^p4p\FS}. (64) 

<2,P,P' 

The constraint from the first matrix element tells us 
that the intermediate state can differ from the Umklapp 
state |m = 2, N) in only two ways since the density op- 
erator can either create a new particle-hole pair, or move 



(a) (b) 




FIG. 4: The states (a)-(d) indicated above are the only in- 
termediate states that give non-zero contribution to the first 
order in V in the perturbative expansion of the form factor in 
Eq. (|63p . This form factor is used to determine A2 to leading 
oder in V. Contributions to the form factor from the states 
(a) and (b) pictured above are calculated in Eq. (|65p and the 
contributions from states (c) and (d) in Eq. (|67p . 



an existing hole or particle to a different 'spot' in mo- 
mentum space. 

At this stage there are two possible types of interme- 
diate states I a). One possibility is a state with only one 
particle-hole pair on top of the ground state such that 
this pair is one of the ones contained in the configuration 
\m ~ 2,N). In this case the density operator creates 
the second excitation of roughly p « 2k p- The second 
possible type is an intermediate state which contains two 
particle-hole pairs. Such a state can at most differ by one 
particle or hole from the configuration |m = 2, N) since 
we can use the density operator to move this particle or 
hole to give the correct final configuration. 
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In the case where \a) contains only a single particle- 
hole pair, the only way the second matrix element can 
be non zero is if the action of the V creates one parti- 
cle hole pair of momentum p ^ 2kp and the other of 
momentum 0. This cannot be done since the operator 
is explicitly momentum conserving. Thus we must only 
consider intermediate states with two particle hole pairs 
where one particle or hole may differ from the configura- 
tion of |m = 2, A^). 

We can further narrow down the intermediate state \a) 
to the two shown in Fig. IDJa), (b). The reasoning is the 
following; let us try to generate a state that differs by 
one particle from the desired Umklapp state. To do so 
we use one of the creation-annihilation pairs to gener- 
ate a hole at either —kp or —kp + 2-k / L and a particle 
at either kp + 2tt / L or kp + Att/L. The second annihi- 
lation operator must then make the hole that was not 
made by the first. This choice along with the momen- 
tum conserving constraint fixes the last particle to be at 
approximately —3k p. The above procedure will generate 



eight terms - for each of the states in Fig.|4l^a), (b), there 
are two choices for where one of the creation operators 
can act and two independent choices for where one of the 
annihilation operators can act. The second creation and 
annihilation operators are fixed. This corresponds to a 
total of four ways to connect the ground state to each of 
the two states, giving a total of eight terms. 

If on the other hand we try to generate a state that 
differs by one hole from the Umklapp state, we find that 
there is no way to do so given the kinematic constraints 
imposed by the interaction term, i.e. using a similar ar- 
gument to the case of the deviant particle, we find for the 
state differing by one hole the only way to conserve mo- 
mentum is to have the annihilation operator act outside 
the Fermi sea. Thus we exhaust all the possible interme- 
diate states that give non zero contributions for both the 
matrix elements. 

We may collect the various contributions to Eq. ([M)) 
for the states shown in Fig. 2] (a) and (b) : 



-kj, {-kp + 2TT/L)^ {kp + 2n/L)^ {-ikpf ^kp f tt 
Efs — — 777T ^ KIT? 7777 77177 — = 177~ \ ' 



2M ' 2M 2M 2M M ' L 



4fcF / 
' M \ 

M{V{2kp + ^)-V{2kp)) M{V{-2kp-^)-V{-2kp)) 
" 8L^kp{kp + ^) 8L^kp{kp + f ) ' 

i-kp)^ {~kp + 2t: /Lf {kp + Air/Lf {~3kp ~ 2tt / L)^ _ Akp f Stt 27r 



2M 2M 2M 2M M \ L kpL^ 

M{V{2kp + ^) - V{2kp + ^)) M{V{-2kp - ^) - V{-2kF - ^)) 



L 

— kTL^ I oi.-r.p ,r.p -r — -r 



= ^ ^ / ^ , ^ V^^- (65) 

SL^kp kp + ^ + M^] SL^kpl'--'^^^^ 



Next we consider the second term that is linear order in 
V(S): 

S 2L^(S„L-E.) '°lS:/t..-/v|FS) 

|a> q',k' 

(m = 2,N\J2 V{-qHl+,i'l,_^^Pp4p\a). (66) 

We find that for the first matrix element to be non 
zero, I a) must contain only one particle hole excitation 



over the ground state which the density operator will 
remove. Furthermore, because V conserves momentum, 
I a) must have momentum Akp. The interaction term 
must act on the state \a) in the following way; one of 
the pairs of creation-annihilation operators must create 
one of the paired particle-hole excitations of momentum 
« 2kp as found in |m — 2,N). The other is automati- 
cally constrained to change the momentum of a particle 
by « 2k p. The combined action of the four fermion op- 
erators on I a) must generate |m = 2,N) in order for the 
matrix element to be non zero. The allowed interme- 
diate states that give a non zero contribution from the 
term are iUustrated in Fig. |3] (c), (d), and give the 
following contributions: 
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Em=2,N — E^c 
Em=2,N — Eld 



4d 



{kF + 2TT/Lf (fcF + 47r/L)2 {-kpf {-kF + 2Tr/L)^ {-kpf (3/cf + 47r/L)2 



2M 2M 2M 2M 2M 2M 



M \ ^ L 

M{V{2kF + x) - yC^kp)) M{V{-2kF - - V{-2kF)) 

SL^kFikF + i) 8L2fcF(fcF + f) ' 

(fcF + 27r/L)2 (fcF + 47r/L)2 (-fef)^ (-fcF + 27r/L)2 {-kF + 2Ti/Lf [ikF + Qir/Lf 



2M 2M 2M 2M 2M 2M 



AkF ( 37r 27r 
+ — ' 



M \ L kpL 

M{V{2kF + x) - ^(^fcf + x)) ^^(^(-2fcF - + - x)) 

8L2fc^(fcF + ^ + ^) " ^L^kF{kF+^-^ + ^) 



(67) 



Collecting these terms gives the following result to lead- 
ing order in 1/ L,V 



(m = 2,7V|/5(0)|G5) 



2MiT^ 



[V'{2kF)-kFV"{2kF)]. 



This means that for A2 we have 
2{poL) 



(27r)8po 



|(m = 2,iV|p(0)|G5)|^ 



[r(2fcF)-A:FV^"(2A;Fr- (68) 



F. Calculation of prefactor Ao,-(fc) 



We can identify the prefactor Ao-{k) using: 



|(fc,7V-l|V(0)|7V)|2 = l + (l_^,__)log 



2-K / .(■}\ /_ I9A . . ( 2tx 



f (1 + <^ + 0(y^)) (1 + (1 - M^^) ) log (^^ j + 0(log^ F 
(1 + 4^1 + (1 - M^^l) log (^) + F^)) , 



(2) _ , [y(0) - V-(fcF + fc)]^ rv? [y(0) - ^(fcF - fc)]^ 

4^2 (fcF + fc)2 47r2 (fcF-fc)2 ' ^"""^ 



where in the second line of the above expression the first 
pair of parentheses contains the perturbative expansion 
of the prefactor and the second contains the expansion of 
the power law. The superscript (2) on a term indicates 
that the term is of 0(V^2) accuracy. 

Similarly to Eq. (l28l) , the form factor required to obtain 
AQ^_{k) is 

(fc, N - \\mm = ^ - ii E V 1^)- (70) 



The state \k,'N — 1) contains a single hole of momen- 
tum k. The perturbative calculation of the form factor 
in Eq. (1701) is complicated by the fact that at the zeroth 
order, i.e. for the non-interacting Fermi gas, the form 
factor is already non-zero. Consequently any dependence 
on the pair potential V enters as a term beyond leading 
order. This is problematic because terms that arise af- 
ter the first non-vanishing order in the expansion of a 
generic parent form factor make two types of contribu- 
tions. They not only contribute sub-leading corrections 
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to the prefactor, which we require, but also generate pow- because there are no 0{V) contributions to the form fac- 
ers of log(L) which give rise to the non-trivial power law tor. 
in L. 

Let us first expand the ket vectors in Eq. ([701) up to 
second order in V and then show that this is required 

I 



77 J^k,N-l — -CJa (-t^k,N-l — -l^a)(J^k,N-l — -t^fS) 

I") I"), 1/3) 

^ (fc, iV - 1| Wy|fc, N - 1) (°)(a|y|fc, N - 1) W , ^ K^^IMV-^W 

|q) ' |q) 

^ |(a|y|FS)p 
I") 



Clearly there is a trivial zeroth order contribution, 

-^(fc,7V-l|(0V3fc|FS) = (72) 

where the annihilation operator removes the particle of 
momentum k from the ground state to connect it to 
|fc,iV- 1>("). 

Furthermore there can be no 0(V) terms. There are 
two sources of first order terms corresponding to keep- 
ing a first order term in one of the ket vector expansions 
and the zeroth order term in the other. Moreover, due to 
the momentum conservation constraint imposed by V the 
intermediate states corresponding to these two groups of 
terms will have a total momentum of either (same as 
|FS)) or —k (same as \k,N — l)^*^-*). Finally, these in- 
termediate states must be connected to \k,N — 1)^°^ or 
|FS) respectively, by the action of a single annihilation 
operator. The only zero momentum state that can be 
connected to |fc,iV- 1)'°^ by a single annihilation opera- 
tor is |FS) since the annihilation operator has to remove 
momentum |fc| < kp i>y creating a single hole, necessarily 
carrying momentum k. On the other hand the same ar- 



gument holds for intermediate states of momentum —k 
connecting the ground state by the action of on the 
left, i.e. the only state that can be connected to |FS) in 
this way is \k,N — 1)'"^. Since the sum over intermedi- 
ate states specifically forbids them from being identical 
to the unperturbed state, we obtain no contribution of 
0{V) to the form factor. 

There are three ways of generating 0(V^^) terms: we 
may either pick 0{V) terms in the ket vector expansions 
of both states in Eq. ((TO)) , or pick an 0{V'^) term from one 
ket vector and an unperturbed term from the other (there 
are two ways of doing this) . However not every term will 
give a contribution. Some terms can be discarded for the 
same reason the 0{V) terms drop out. This is because 
the structure of the matrix elements appearing in these 
terms is the same as in the 0{V) terms - i.e. they require 
the matrix element of the annihilation operator between 
either |FS) or \k,N — 1)''^^ and an intermediate state to 
be zero with the constraint that the intermediate state 
must have the same momentum as, but cannot be equal 
to either of those states. Thus we need only consider the 
following terms: 



15 



p _ l">:l^> P" Pl,Pi,«l P2,P2''?2 

(73) 

8L5/2(iJps _ £;„)2 ' ^^^^ 

EE E (fc,^-l|'"^V|FS)|(fc,A^-l|(°V(gi)^^^+,^V3^,_,,^p;V'pJa)P 

where 

(/c,iV-l|(")4|FS) - -^+i^i+F2+i^3 + 0(y4). (76) 



Let us now consider each of the terms in Eqs. ([75)1 
- ([75]) on the basis of allowed intermediate states. We 
observe that the only admissible states |a) in Eq. ([73)) 
are ones containing a net zero momentum pair of particle- 
hole pairs and a hole of momentum k. This is because the 
third matrix element in Eq. ([75)) will lead to states 
with two net-zero momentum particle-hole pairs, while 
the second matrix element can indeed admit two kinds of 
states. However, the matrix element of the annihilation 
operator will only connect |a) and that have identical 
pairs of net zero momentum particle-hole pairs, since the 
annihilation operator can then create the additional hole 
at k in (3. Thus we may write the total contribution due 
to Fx as 



Fi = 

E E E 



2[V{q)-V{p-p' + q)f 
4L5/2(A£;)2 



AE = -^{q+{p^p')). (77) 
m 

In Eq. ([77)) we enumerate all possible ways to gen- 
erate an intermediate state with two net zero momen- 
tum particle-hole pairs. We can understand the limits 
on the sums in the following way: the first hole can be 
placed anywhere in the filled Fermi sea (except at the pre- 
existing hole of momentum k). With this hole in place, 
the momentum transfer to the corresponding particle (to- 
tal momentum oip + q) needs to be large enough that the 
particle is created "outside" the filled Fermi sea. We first 
consider a positive momentum transfer and have that the 
transfer q must be larger than kp—p. This sets the lower 
limit on the second sum. Finally we consider the second 
particle-hole pair. For the hole of momentum p' to lie 
within the Fermi sea and simultaneously have its corre- 
sponding particle of momentum p' — q outside the sea, we 
note that p' can have at most a momentum of q — kp pro- 



vided that q — kp < kp. Ifg exceeds this latter condition 
then p' may lie anywhere in the sea as long as it is not 
equal to p or fc. If we then consider how the intermediate 
state corresponding to some allowed choice of p,p' ,q can 
be created with a negative momentum transfer we find 
that setting q ^ p' — p — q also generates the same state. 
Moreover we may also relabel p ^ p',p' ^ p and repeat 
the same argument above, leading to an additional fac- 
tor of 2 multiplying the sums. The relative signs of these 
terms is fixed by the canonical anti-commutation rela- 
tions of the fermionic operators. Furthermore the only 
way to connect states a and P in Eq. ([75)) with the an- 
nihilation operator matrix product is if the particle-hole 
pairs in the two states are identical, leading to the ap- 
pearance of the square in Eq. ([77)) . 

A very similar argument to the one above allows us to 
write down an expression for F2. The key difference is 
there is no restriction disallowing the holes from having 
a momentum k since the intermediate states are not in- 
volved in the annihilation operator matrix element. Thus 
we have 



Fo 



kp oo min[feF,i3— fcp] 

EE E 

P^ — ^F q>kF—p p'/p, 
p' > — fcp 



2[V{q)~V{p-p' + q)f 



AE = -^{q + {p^p')). (78) 
m 

Lastly, for i^a we have two groups of terms. On the 
one hand akin to Eq. ()77p we may have intermediate 
states with two particle-hole pairs. Here too there is a 
restriction on where the holes can be since there is a pre- 
existing hole of momentum k. On the other hand, we 
may also have states where the hole of momentum k has 
been 'moved', followed by the creation of an additional 
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particle-hole pair. 



kF 



oo Ta\Ti[kF,q—kF] 



p^k, qykp-p p'^k,p, 

P>~kF p'>-kF 

p<k —kp—p+k 



2[V{q)-V{p-p' + qf 
8L5/2(A£:)2 



E E 

p—~kF p' — — kF 
kF kF 

E E 

pyk p' >kF — 'p-\-k 



2[V{k^p)-V{k-p')f 
8i5/2(A£;)2 

2[V{k-p) -Vjk-p')]'^ 
8L5/2(A£;)2 ' 



AE = -^iq+ip-p')), 
m 

/^E^-{p + p' -k)-?^. 
m m 




kp + k — p 



ki-A) 



P 



FIG. 5: The various regions of integration corresponding to 
summing up the contribution to the form factor from states 
that have two holes and one particle excitation. 



(79) 



Thus upon summing these contributions we are left 
with the following terms: 



E^^ - 



(1) _ 



(T40X + T^cl + ^3c2)i 

p<k —kp—p+k 



rr.(2) _ 

^3cx — 

AE = 

AE = 



E E 

p— — kF p' — ~kF 
kF kF 

E E 

p>k p'ykp—p+k 



00 min[q—kF.kF] 

E E 

q>kF-k p'^k, 
'p'~> — kF 

-±{q + {k-p')), 
m 

— {p + p -k) 

771 771 



[V{k~p) - V{k-p')f 

4L5/2(A^)2 

[V{k-p) - V{k~p')Y 



4L5/2(Ai;)2 

[V{q)-V{k~p' + q)f 
4i5/2(A£;)2 



(80) 



The notation above is suggestive of the fact that we are 
considering contributions to the form factor due to inter- 
mediate states with 4 excitations (2 particle-hole pairs) 
and 3 excitations (a particle-hole pair and an unpaired 
hole) separately. 

Let us consider first: 



2 p<k p' <-kF + k~p r , , , ,>,2 



3cx 



4L5/2 



p——kF p' — — kp 



{k-pY{p' -k) 



fe— A — feF+fc— p 

-2 dp f dp' [V{k-p)-V{k-p')f 



2L1/2 7 27r 7 27r {k~pf{p'-kf 

-kp+k —kF 
2 

^2 p' = -..+fc-p-2./L j^^^ _ _ _ ^,^j2 



2L5/2 



E E 

p—k—A p' — — kp 

2 2 

dp 
2^ 



(fc-p)2(p' -fc)2 

dp' [V{k-p)^V{k~p')]^ 



4L1/2 7 27r 7 27r {k-pY{p'-kY 

—kF —kF 

Starting from the initial sum, we have separated the 
contributions from three different regions over which we 
need to sum the right side of the above equation. We 
have taken the continuum limit to express the sums as 
integrals in these regions, however, we have separated the 
contributions in regions (1) and (2) when one of the holes 
is within the region of momentum space, (fc— A, k) (corre- 
sponding to the lightly shaded triangles in Fig. ([5])), and 
maintained these contributions as discrete sums. This is 
because the energy denominator will become badly be- 
haved in this region and upon integration yield a loga- 
rithm of L; however we require more than logarithmic 
accuracy since such corrections may enter the prefactor. 

While our expressions are formally correct for a fixed 
A. We wish to demonstrate that the limit as A — > 
of the above expression exists and subsequently obtain 
an analytic expression for it. In order to do this we will 
approximate the integral when the range of p nears A: — A 
and obtain the leading order A dependent contribution. 
We will then carry out the discrete sums and obtain the 
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leading order A dependent contribution as well as any 
non-vanishing constants. We will show the existence of 
the A limit by demonstrating that all divergent A 



J 



dependent contributions vanish. 

After approximating the behavior of the divergent in- 
tegral in Eq. (|8T|) at the upper bound we are left with 



167r2fc2Li/2 



dx 



dy 



[V{k{l - x)) - V{k{l - y))]^ m,^[V{0) - V{k + kp)]^ 



-k p / k 



-kp/k 



(l-.T)2(y-l)2 



(log(A)) , (81) 



r 



where the notation and meaning of the special principal 
value integration denoted by V+ is defined in the ap- 
pendix, see Eqs. ((CT|) - ((C2]) . 



J 



Meanwhile we may evaluate the sum in Eq. (jSip ex- 
plicitly using the identities and expansions of the poly- 
gamma functions, Eqs (|Bll) - (jB3l) : 



p—k-—27r I L p — — kp--\-k—p—27T/L 



E 



E 

'p'=-kF 



'M2 



p=k-A 

m2[T/(0) - y(fc + fcF)]2£3/2 

327r4((iV-l)/2 + nfe)2 
m^[V{Q)~V{k + kF)]'^ 



[V{k-p)-V{k-p ')] 
{k-p)'^{p' -kf 



log(A) - log ( ^ 



- IE 



- IE 



(82) 



r 



where 7^ k, 0.5772 is the Euler-Mascheroni constant. 

We note how the logarithmic divergence appears with 
opposite sign in the finite sum and the boundary of the 



J 



integral and thus drops out of the final answer. We there- 
fore obtain an expression for the term 



(1) 



3ex 



7^4 



1 —kp /k+l~x 

dx f dy [V{k{l - x)) - V{k{l - y))]2 



4L1/2- ^ J 2tt J 27r 

l-kp/k -^kp/k 



fc2(l -2;)2(y- 1)2 

dp' [V(k - p) - V{k - p')]2 m2[y(0) - V(k + kp)]^ 



— fc^-|-fc —k^-{-k 

9 ^ ^ 
m f dp f 

4L1/2 / 2^ J 2tt {k-pf{p'-kf 



-k. 



87r2Li/2(fc + fcF)2 



log(^ 



1e]. (83) 



r 



A similar calculation yields for the second term sum- excitations 
ming contributions from intermediate states with three 
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T. 



(2) _ ™' r dx f dy[V{k{x-l))~V{k{y-l))]^ 



3cx 4^1/2 - J 2tt J 2tt fc2(l-a;)2(y-l)2 





dx 


- / 






1 






kp 


f dp 


I 


J 27r 




2 


kp- + k 

2 



' ^og f ^ ) _ 7s ) . (84) 



{k-py{p'^ky 87r2Li/2(fc^ _ fc)2 V 

I 

The sums from the contribution of intermediate states haved and can be immediately interpreted as integrals: 
with two pairs of particle-hole excitations are well be- 



m2 ^ "'"'^'!^'-'''^ [V{q)-Vik-p' + q)]^ 



4L-V2 



5/2 X! X! 



qykp—k p'>-~kf 



q^k-p' + qy 



m 



167r2Li/2 



/ dp 

kp—k J —kp 



q^{q + k ~ p')2 
I 



, r"^, ,[Viq)~Vik~p' + q)]^ 
aq / dp 

2kF <J —kp 



q^{q + k~p'Y 



(85) 



Thus we obtain the full contribution to the form factor 
up to order V"^: 



{k,N ~l\i;mN) 



1 



2L1/2 



-fcp / k-\-l—x 



dx J dy 

l-kp/k -kp/k 
2 

kp /k 



[V{k{l-x))^V{k{l~y))] 



dx 



dy 



/c2(l-x)2(j/- 1)2 

[V{k{x-l))-Vik{y-l))r 



1 kr/k+l-x 

~-kp -\-k — k p- -\-k 



fc2(l-a;)2(y- 1)2 

kf kp 



/M2 



dp I dp 



2kp 



'kj- 



, [ V{k-p)-V{k~p ')] 
2{k-pY{p' -kf 



dp / dp' 



F + k kp +k 



, [ v{p -k)~ v{p' - k)Y 

2{k-pY{p' -kf 



Ikp-k^J-kp " 2q^q + k-pr 
/ [V{0)-V(kF + k)]^ [VjO) ^V{kF-k)]^ 
^[ {kp + k)^ ^ {kp-k)^ 



IE 



4^2 



[V{{))-V{kF + k)Y [V{Q)-V{kF-k)f 



{kF + kf 



{kF - kf 



. 2q^{q + k-p'Y 



log(^ 



(86) 



Comparing Eq. (1861) with the field theory prediction, 
we obtain the leading correction to the prefactor Aq.- {k)'- 



19 



dx 



dy 



1-kp /k 



-kp/k 



[V{k{l^x))-V{k{l-y))Y 
fc2(l-a;)2(y- 1)2 



l+kp/k 



kp I k 



+V- J dx 

1 kp/k+l-x 

— kf+k —kp-\-k 



dy 



[V{k{x-l))-V{k{y-l))f 
F(l-a;)2(y-l)2 



dp I d/''t7'^:J^'-'?' + / dp 7 dp' [^(^ - ^) - - 



2(A:-p)2(p' -fc)2 



2{k-pY{p' -kf 



2kf 



fcp-fc 2q^{q + k — p'Y jjfcr ^'-fcF 

[y(o) - + fc)]2 [y(o) - y(/cF - k)f 



2q^{q + k-p'f 



{kp + ky 



{kp - kf 



IE 



(87) 



IV. SUMMARY 

We present below a few results for the nonuniversal 
prefactors of correlation functions and dynamic response 
functions (see Eqs. dH)-© and Eq. The results 

of the perturbative calculations described in the previ- 
ous section produce the following leading order results 
for the prefactors (see Eqs. (EH) , ([Ml) , (El), ©, and © 
respectively): 



A2 



A'P{V(kF + k)- Vjkp - fc))^ 
87r3(fc2,-fc2)2 ' 

AP{V{2kF)-V{kF-k))^ 
327r3fc2,(fc-/ci.)2 ' 

M^Vik-kp)^ 
Sn^ik-kp)^ ' 
M^PoV'{2kp)^ 

M^poiV'i2kp) ~ kpV"{2kp) f 



327r 



10 



(89) 
(90) 
(91) 
(92) 



We see that all the prefactors above vanish when V{k) 
is constant, satisfying the expectation that a contact 
interaction V{r) oc 5{r) between the spinless fermions 
should not affect any observables. We were also able to 
reproduce the known values of the prefactors Ci,A2 ^ 
for the Calogero-Sutherland model, with our perturba- 
tive fermionic results, Eqs. (22), (23). We note that the 
results above would be hard to obtain using the conven- 
tional infinite size diagrammatic technique due to log- 
arithmic divergences [Tl|, |l3], but the consideration of 
finite size scaling allows us to obtain these answers using 
essentially only undergraduate quantum mechanics tools. 
To lowest non- vanishing order, we see that the prefactors 
Ao.+(^): ^i,±(fc) are oc . In general, to lowest non- 



vanishing order ^„>i_±(fc) will be oc V'^'^ , although the 
complexity of the expressions grows. To obtain correc- 
tions beyond leading order to Ari,±(fc) (and the first non- 
trivial correction to Ao._(A:),Ai and Co), one needs to 
carefully separate terms in the perturbation series that 
diverge as powers of log(L) from the relevant correction 
to the prefactors, which we have demonstrated in Sec- 
tion IIIF above. 

For the Lieb-Liniger model, we obtain analytic ex- 
pressions for the prefactors of correlation functions non- 
perturbativel y. We relate them to the exactly known 
form factors |3l| which have been used recently for nu- 
merical evaluation of dynamic response functions 36]. 
In Fig. [6] we show results for a few prefactors for differ- 
ent interaction strengths (details will be presented else- 
where [13]), and some known limits are also plotted for 
comparison. As a further check we have also successfully 
reproduced the perturbative fermionic result for the pref- 
actor A2 [see Eq. ([Ml) ] of the Cheon-Shigehara model (38i |. 
which is dual to the Lieb-Liniger model and has the same 
density correlations. It should be noted that recently 
the correlation functions for the Lieb-Liniger model have 
been treated without field theoretical considerations but 
rather by directly wor king with the microscopic details 
of the theory (see Ref. [331), and the results of this work 
coincide with ours. 

To summarize, we developed a general approach to cal- 
culating nonuniversal prefactors in static and dynamic 
correlation functions of ID quantum liquids, by relating 
them to the finite-size scaling of the matrix elements of 
the corresponding operators. To find a given prefactor, 
only a single matrix element (form factor) between the 
lowest energy states needs to be evaluated, see Eqs. pT|) - 
P^ . (pS)) . Moreover, the method does not rely on the 
integrability of a model. To demonstrate our approach, 
we calculated some prefactors in static and dynamic cor- 
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' \ Tracy- Vaidya Limit (0.5214) 



FIG. 6: (Color online) Non-perturbative results for the Lieb- 
Lmiger model of ID bosons: 27r^yli (dashed black), 5127r^°A2 
(solid green) and Bo (dot-dashed blue), —Sii^Bi (dotted red) 
as functions of the Luttinger liquid parameter K. In the limit 
of strong interaction {K — >■ 1) our expressions for Bo and Bi 
agree with the known values 7], while Ai — >■ 1/2tv'^,A2 — >■ 
are in accordance with the density correlator of the free 
Fermi gas. We also match Bq in the weakly interacting regime 
{K ^ 1) to Popov's result (dashed line) [^, and show some 
numerical results (crosses) [40| ]. 



Spectral Weight /^^ 



A 



/ r(^R+3) , 



p = 0,g = -2 p=l,5 = -l 



AEocl/L 

FIG. 7: (Color online) Schematic showing how spectral 
weights are distributed amongst multiplet of states which are 
degenerate within oc 1/L accuracy. In the inset, we illustrate 
how spectral weights are distributed within a multiplet for 
n+ = 2, once the degeneracy is lifted by oc l/L^ corrections. 



where 



C{nr,HF..B) = 



TOiF,R)^{n+ + 1) - 



(A2) 



and f{{pi,qi}) is defined in Eq. (pO|) . The rules for n_ 
are analogous. 



relation functions for weakly interacting spinless fermions 
with an arbitrary pair interaction potential, see Eqs. ([88])- 
(192]). 

Methods of Section II applied to an integrable model, 
allow one to obtain exact expressions for prefactors of 
the response functions, at arbitrary interaction strength. 
The results of the perturbative calculation detailed in 
Section III agree, in the proper limit, with the results [s^l 
obtained by the methods of Section II for the exactly 
solvable Lieb-Liniger model, see Fig. [5] 
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Appendix A: Multiplet Summation Rule 

Let us now provide an illustration of the multiplet sum- 
mation rule, Eq. ((22)) of the article, for a few simple cases. 
The sum rule is given by 

E \faP^}■A<l^})\^ = C{n+,^iF^J,), (Al) 



1. n+ = 1 

For = 1, the only state which contributes to the 
multiplet is p — 0,q = —1. Thus the left hand side of 
Eq. (|A1I) is given by 



|/({0,-l})p = |/+(0)x/-(-l)p 



{-^/Wm - I)--. 
/ (^//^- 1)(V7^- 2). 



(A3) 



The right hand side of Eq. (|Aip is given by 
G(n+ — \-,fiFM) 



[PF,R - 1) 

l/({o,-i})r 



2. n+ = 2 



^'■F,R 

(A4) 



For n+ = 2, there are two states: p = l,q = — 1, 
and p = 0,q — —2. Consequently the left hand side of 
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Eq. (jAll) is given by 



|/(+l; -1)1^ + 1/(0, -2)1^ = 

2 



i/+(i)r(-i) 



i/+(o)r(-2) 



fJ-F,RipF,R + 1) 



(A5) 



And the right hand side of Eq. (jAl|) is given by 

'{t^F,R + i){t^F,R){t^F,R - !)• 



IJ-f,r{iJ-fm + 1) 



r(3)(AiF,7? - l)... 
=|/(+l;-l)P+|/(0,-2)p. (A6) 



3. n+ = 3 

There are three states in the niuhiplct ?i+ = 3: p = 
2,(7 = -l,p = 1,9 = -2 and p = 0,q = -3. The left 
hand side of Eq. (IA1[) is given by 

|/(+2;-l)|2 + |/(+l;-2)|2 + |/(0;-3)|2 
1 



3d 



r(MK + 3) 



6r(Mfl,) 

Thus once again the summation rule holds. 



(A7) 



Appendix B: Identities involving Poly-Gamma 
functions 



We have a few explicit values of the polygamma func- 
tions: 



^(o)(l) = -^s, 



V.«(l) 



6 



(B2) 



where ~ 0.5772... is the Euler-Mascheroni constant. 
We also have the following asymptotic expansions: 



2z 



hm ^(^\z)^l + -L + o(^ 



z 2z2 



(B3) 



Appendix C: Principal Value Integrals: Interior and 
Edge Singularities 

We define three distinct types of principal value inte- 
grals as follows. 



X — C 5-i.O 



X — C 



X — c 



P t dxl^ = hm (t 1^+ J [a] log(5) I , 
I'dxIH = lim ( [' ' M,^ f{h)\og{5)\ , 



X 



- b s^o 



X 



p,i dx-im^^p.r^-p 



{x-a){x-b) '^Ja x-b x~a 

(CI) 



Note also that when we have a multiple integral where 
all but one of the integrals is analytic and do not need to 
be evaluated in the special sense above, we will use the 
notation 



The poly-gamma function is defined as (z) — 
log(r(a;))|a;=2. We have the following summation 
identity for ?/''^"^(z): 



V i 

^ (fc + a)"+i 



(^(")(Ar + a) _^W(a)).(Bl) 



a -/ c 



X — 6(a) 



(C2) 



where we mean that we can freely carry out the inner in- 
tegration, but the final integration needs to be performed 
as in Eq. JCl]). 
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